Local exclusion and Lieb-Thirring inequalities 
for intermediate and fractional statistics 



Douglas Lundholm 1 '* and Jan Philip Solovej 2 '^ 

1 Institut Mittag-Leffler, Auravagen 17, SE-182 60 Djursholm, Sweden 
2 Department of Mathematical Sciences, University of Copenhagen 
Universitetsparken 5, DK-2100 Copenhagen 0, Denmark 



Abstract 

In one and two spatial dimensions there is a logical possibility for 
identical quantum particles different from bosons and fermions, obey- 
ing intermediate or fractional (anyon) statistics. We consider appli- 
cations of a recent Lieb-Thirring inequality for anyons in two dimen- 
sions, and derive new Lieb-Thirring inequalities for intermediate statis- 
tics in one dimension with implications for models of Lieb-Liniger and 
Calogero-Sutherland type. These inequalities follow from a local form 
of the exclusion principle valid for such generalized exchange statistics. 

MSC2010: 81Q10, 81S05, 35P15, 46N50 

1 Introduction 

Fundamental to many-body quantum mechanics is the notion of identical 
particles and associated particle statistics. Many of the remarkable quantum 
phenomena observed in nature are direct consequences of the fact that iden- 
tical particles in three spatial dimensions are either bosons, obeying Bose- 
Einstein statistics, or fermions, obeying Fermi-Dirac statistics. For example, 
the amplification of light in lasers by superposition of photons is possible 
because they are bosons, while the periodic table of elements arises from 
the sequential filling of atomic shells by electrons which are fermions. This 
binary classification has been known since the early days of quantum theory, 
and follows logically from the requirement that an iV-body wave function 
^(a5i,...,jCjv), modeled as a square-integrable complex- valued functioijll on 
the iV-particle configuration space (R 3 )^, or alternatively as an element 
of the tensor product of one-particle spaces (g) w L 2 (]R 3 ), has to be either 



*Work partly done while visiting FIM, ETH Zurich, (e-mail: lundholm@math.ku.dk) 
^Work partly done while visiting Institut Mittag-Leffler. (e-mail: solovej@math.ku.dk) 
1 We will for simplicity always model particles as non-relativistic and scalar, i.e. with- 
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symmetric (corresponding to bosons) or antisymmetric (corresponding to 
fermions) under the action of permutations of the particle labels — which, 
if the particles are indistinguishable, must leave the observable probability 
\ip\ 2 invariant. However, it was discovered in the 1970-80's [5^ 1 |25 | [19 ], I5S] 
that this simple picture for identical particles and statistics is logically in- 
complete and should be replaced by a more general framework, which in 
three (or higher) space dimensions nicely reduces to the above binary pic- 
ture, but in the one and two dimensional cases actually admits a wider range 
of possibilities. Identical particles classified via this generalized framework, 
but falling outside the usual boson/fermion classification, were said to obey 
intermediate or fractional statistics. Although elementary identical par- 
ticles live in three space dimensions and hence must be either bosons or 
fermions, these more exotic one- and two-dimensional possibilities have re- 
cently become more than just a mathematical curiosity, with the advent of 
trapped bosonic condensates [M1H2] and quantum Hall physics [2l[37]> and 
thereby the discoveries of effective models of particles (or guosi-particles) 
which seem to obey these generalized rules for identical particles and statis- 
tics (see [SI [13 E31 ESI SH S3 SSI EZ] for reviews) . 

In two space dimensions the usual requirement, that the phase change 
of the wave function which arises from an interchange of two identical par- 
ticles (or, in the more precise general framework, from a continuous simple 
interchange of two particles) needs to be either +1 (boson) or -1 (fermion), 
is no longer valid. This interchange phase is instead allowed to be any fixed 
unit complex number e m7r . The corresponding particles are called anyons, 
and are alternatively parameterized by the real-valued statistics parameter 
a G (—1, 1]. A standard way to model such anyons is by means of bosons in 
M. 2 together with magnetic potentials Aj(x) of Aharonov-Bohm-type with 
strength a between every pair of particles (see (p7|) below), giving rise to the 
correct phase of the wave function as the particles interchange or encircle 
each other |25l I58j . The free kinetic energy operator for N anyons is thus 
given by (see (|27|) ; we use the unit conventions H = m = 1) 

1 - 

3=1 

and acting on completely symmetric wave functions tp 6 (££)^ m L 2 (IR 2 ). 

In one dimension, the above geometric picture of continuous interchanges 
of particles breaks down — particles have to collide in order to be inter- 
changed. In a quantum mechanical context this necessitates the prescrip- 
tion of boundary conditions for the wave function at the collision points 
r := Xj + i — Xj — > + (i.e. at the boundary of the proper configuration 
space X := {x 6 M. N : x\ < X2 < ■ ■ ■ < x^}, in this case with codimension 
one). There are two ways in which one may arrive naturally at a set of 
boundary conditions. One approach |25[ H], which is generally referred to 
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as Schrddinger quantization, starts out with the proper configuration space 
X and imposes a general (Robin) boundary condition of the form d r ip = rjip 
at r = + , with r/ £ M an arbitrary but fixed statistics parameter. The 
case 77 = corresponds to bosons (Neumann b.c.) and r] = +00 to fermions 
(Dirichlet b.c). Such boundary conditions can, if ip is extended from X 
to M. N and modeled symmetric about r = 0, alternatively be encoded by a 
delta potential 2r]5{r), and the resulting model for ID intermediate statis- 
tics, represented with a modified kinetic energy operator 



and acting on completely symmetric wave functions ip 6 (^)^ m L 2 (R) , is 
equivalent to the Lieb-Liniger model [30] for a gas of bosons on M. with delta 
interactions. 

Another, inequivalent but just as well-motivated, approach to ID inter- 
mediate statistics is to consider Heisenberg quantization [261 EH Ell EZ] in 
the sense that one starts out with the classical phase space M N x VL N and 
its observables (both should be symmetrized w.r.t. relabeling of particles) 
and considers representations of a corresponding algebra of operators. One 
finds, in the two-particle case, that the relevant representations are labeled 
by a continuous parameter a £ K, and that in the standard coordinate 
representation this amounts to imposing a boundary condition ip(r) ~ r a 
for the wave function as the relative coordinate r — > + (this relation with 
boundary conditions is more complicated and will be elaborated in Section 
12.11 below, but the restriction a £ (—1/2, 00) is always required for square- 
integrability) . In this context bosons correspond to a = and fermions to 
a = 1. Again, using symmetric tp the general boundary condition or choice 
of representation can equivalently be encoded by an interaction potential, of 
the form a(a — l)/r 2 . The resulting model can be extended to the A r -particle 
case and is equivalent to the homogeneous part of the Calogero-Sutherland 
model [HIES], represented with a modified kinetic energy operator 



and again acting on completely symmetric wave functions ip G ® sym ^ (^) • 
One may ask what properties these different types of identical particles 
possess. The fundamental characteristic property of fermions is that they 
obey Pauli's exclusion principle, i.e. the fact that the antisymmetry of the 
wave function implies that no two particles can occupy the same single- 
particle state: <p A <p = for any ip 6 L 2 (M. d ). A powerful consequence of 
this is the Lieb-Thirring inequality [MIES] (see also [53]) for the energy of 




j 



l<j<k<N 
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N (spinless) fermions in an external potential V in ~R d , given by (again, with 
the conventions H = m = 1, and H^H = 1) 



N-l 

> - J2Mh)\ > -C d \V-(x)\ 1+ % dx, (1) 
fc=o 

where V± := (V ± |V|)/2 and C<2 a positive constant. Here \k(h) denote 
the negative eigenvalues (ordered by decreasing magnitude and with multi- 
plicity) of the one-particle operator h = — + V(x). The first inequality 
expresses the Pauli exclusion principle while the second concerns the trace 
over the negative spectrum of h and in some sense incorporates the uncer- 
tainty principle. In fact, the Lieb-Thirring inequality is equivalent to the 
kinetic energy inequality 

N 

T *--= \ [ El V ^| 2 dx > C' d [ p(x) 1+ Ux, (2) 

2 jR dN T~f JR d 



with C' d := d(2/C d ) 2/d /{d + 2) 1+2 / d another positive constant, and p the 
one-particle density function associated to tfi: 



N 



P( x ) : =yZ . . ,Xj^,x,x j+1 ,. . . ,x N )\ 2 T\dx k . (3) 

3=1 kfy 

The inequality (J2J can be interpreted as a strong form of the uncertainty 
principle valid for fermions. 

These inequalities need to be weakened in the case of weaker exclusion. 
E.g., in the situation that up to q E N particles can occupy the same one- 
particle state (sometimes referred to as Gentile intermediate statistics \W\ 
117]). the r.h.s. of ([I]) resp. ([2]) are to be multiplied by q resp. q~ 2 l d . 
For bosons there is no restriction on the number of particles in the same 
state; we can e.g. consider ip = ipo ipo ® . . . ® ipo, where ifo is the ground 
state of h. Hence bosons can be accommodated by taking q = N, and 
the inequalities become trivial as N — > oo. However, for intermediate and 
fractional statistics (in the above sense) the picture is more complicated. 
There have been partially successful attempts to relate also such many-body 
quantum states to one-particle states restricted by some exclusion principle 



2 The currently best known value for C' d is a factor {?>/% 2 ) 1 / d smaller than its semi- 
classical value with Cd = f 2 ir) d Jk<*(1 ~ l£| 2 )+^£ BB- ^ ls conjectured that for d = 3 the 
inequality ([2]) for fermions holds with exactly the semiclassical constant and the Thomas- 
Fermi expression on the r.h.s. |34] , 
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[20|. [22l I18j . but so far no general picture has emerged. The difficulty has 
to do with the fact that these generalized exchange statistics are naturally 
modeled using interacting Hamiltonians, hence leaving the much simpler 
realm of single-particle operators and spaces. 

In [38] we found that anyons in M 2 , with statistics parameter a = fi/u 
an odd numerator fraction, satisfy a kinetic energy inequality of the form 



N 

T ^ : =l f yZ\DM 2d x > f p(x) 2 dx, 



(4) 



implying the Lieb-Thirring inequality 

N 



/ (\\D^\ 2 + V( Xj M 2 ) dx > -v 2 C' K [ \V-(x)\ 2 dx, (5) 
Jr 2N V 2 / Jm 2 



for some positive constants C\, C' A (numerical estimates will be given be- 
low). For a general statistics parameter a£l and finite N these inequalities 
hold with 1/v 2 replaced by C 2 N , where 

C a n '■= rnin min \(2p + l)a — 2q\. (6) 

pe{0,l,...,N-2} q& 

However, this expression tends to zero as N — > oo, except when a = n/u 
is an odd numerator (reduced) fraction in which case lim7v->oo C a ,N = 1/v 
(cp. Figure 2 in [38j). This, at first unexpected, difference between odd 
numerator fractions and even numerator and irrational a is further discussed 
in [39], where we also provide some arguments for this difference actually 
being natural. 

In contrast to previous proofs of the standard Lieb-Thirring inequalities 
([1]) and ([2]) , the inequality (jlj) for anyons was proved using a local form of the 
uncertainty principle together with a local form of the exclusion principle, 
valid for anyons with general a and involving the constant C Qj tv in ([6|) mea- 
suring the strength of exclusion. Our inspiration for this local approach to 
Lieb-Thirring inequalities, applicable to interacting systems, was the orig- 
inal proof of stability of ordinary fermionic matter, |12] due to Dyson and 
Lenard (see also [13} 128]). which involved only local, comparatively weak, 
consequences of the Pauli principle for fermions (see Lemma [6] below). 



1.1 Main results 

In this work we consider further consequences of our results in [38] for 
anyons, and extend the above families of Lieb-Thirring-type inequalities 
to identical particles obeying intermediate exchange statistics in one dimen- 
sion, in the above sense. We again start from a local form of the uncertainty 
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principle together with a local exclusion principle of the form (cp. Lemma 
[6] below for fermions) 

J ^x)(r s/H v)(x)dx > (rc-l)j^fp J |^(x)| 2 dx, (7) 

[a,b] n [a,b] n 

where in the Schrodinger case £s is a measure of exclusion depending on r\ 
times the length \a — b\ of a local interval, and in the Heisenberg case £h 
depends only on a (see Figure [Qresp. El below for their exact dependence). 

Theorem 1. Let p be the one-particle density ([3|) associated to a normalized 
completely symmetric wave function ip <G (g)JJ m L 2 (IR) of N > 2 identical 
particles on the real line. 

For the Schrodinger case, with statistics parameter ij > and kinetic 
energy Ts = (ip,Tsip), we have 

T s > C s / ^/~p{x)f P {xfdx, (8) 

JR 

where p is a local approximation of p (which here can be taken to be the asso- 
ciated Hardy- Littlewood maximal function), and Cs is a universal constant 
satisfying 3 • 1(T 5 < C s < 2/3. 

For the Heisenberg case, with statistics parameter a > 1 and kinetic 
energy Th = (ip,Tnip), we have 

T H > Ch^hH 2 [ p(xfdx, (9) 

JR 

for p a local approximation of p given in (|72p (and that here cannot be 
substituted by p or the maximal function), and some universal constant 
1/32 < C H < 2/3. 

In particular, if the density p is confined to an interval of length L, then 
we show that ([9]) implies for the energy per unit length 

T H /L > C H £ H (a) 2 p 3 , (10) 

where p := N/L is the average density of particles. Note that £h(o) ~ ct 
to leading order as a — > oo (see Figure [2]). We can compare this bound 
with Calogero [8] and Sutherland [55] (whose models are exactly solvable 
for certain choices of external potentials), where one finds for the ground 
state energy of a system confined to an interval 

T H /L -« 2 p 3 , (11) 

in the thermodynamic limit N, L — > oo with fixed density p. However, note 
that we do not need to assume a particular confinement potential for our 
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bounds since they concern the kinetic energy Th alone. The possibility 
for considering more general external potentials in the Calogero-Sutherland 
models has also been discussed in the context of Thomas-Fermi theory; see 

[SB EH- 

In the Schrodinger case, the bound ([8]) implies that if p < 7/5 for some 
constant 7 > 0, i.e. if the density is sufficiently homogeneous, then 

T S /L > Cs^v/i^fp 3 . (12) 

We have £s(*) ~ for small t and £s(*) — > as i — >• 00 (see Figure [[]). 
Compare with Lieb-Liniger |30j (again exactly solvable in the absence of an 
external potential), where it is shown that in the thermodynamic limit the 
ground state energy satisfies 

T S /L -> \e{2r } /p)p\ (13) 

2 

with an implicitly defined function e s.t. e(t) ~ t,t <C 1, e(t) — > \,t — > 00. 
We can also compare with [32] , where an energy functional of a form similar 
to the r.h.s. of ([8|) arises in the limit of tubular confinement of a three- 
dimensional bosonic gas. 

The methods we use to prove Theorem [1] are similar to those used for 
anyons in [38J. However, certain new technical complications arise in the one- 
dimensional context, such as a local dependence of the strength of exclusion 
in the Schrodinger case, and the possibility of arbitrarily strong exclusion 
in the Heisenberg case. This is the reason for the more complicated expres- 
sions ([8]) and Q as compared to ([2]) and (JH). There is certainly room for 
improvement in our bounds for the constants Ca/s/h ( we mainly consider 
the forms of the inequalities and their dependence on the statistics parame- 
ters to be of conceptual interest), and for our methods we need to restrict to 
nonnegative statistics potentials, i.e. rj > resp. a > 1. The intermediate 
Heisenberg case a G (0, 1) is certainly interesting (so is a £ (—1/2,0)), but 
presents additional challenges and will not be addressed here. 

In Section [2] we establish the context and notation of the paper, taking 
care to define the relevant operators properly. This is something we have not 
found been discussed in detail in the literature. In particular, we show that 
the formal definition of the kinetic energy form used for anyons in [38] is 
indeed the natural one. In Section[3]we recall the local exclusion principle for 
fermions and anyons, with an application of the latter to an explicit bound 
for the energy of the ideal anyon gas. We then introduce corresponding local 
principles of exclusion for ID intermediate statistics. In SectionH]we deduce 
a local form of the uncertainty principle valid in arbitrary dimensions. These 
local bounds are then applied in Section to prove Theorem [TJ In the final 
section we consider applications of the Lieb-Thirring inequalities dH) and 
([5]) for anyons to the old problem of many anyons confined in a harmonic 
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oscillator potential, as well as the question of thermodynamic stability for 
a system of charged anyons and static particles interacting via 3D Coulomb 
potentials. We end with a discussion on the addition of external potentials 
for the Heisenberg case. 

Acknowledgements. We thank Giovanni Felder, Jiirg Frohlich, Jens Hoppe, 
Edwin Langmann and Robert Seiringer for comments and discussions. Sup- 
port from the Danish Council for Independent Research as well as from In- 
stitut Mittag-Lefher (Djursholm, Sweden) is gratefully acknowledged. D.L. 
would also like to thank IHES, FIM ETH Zurich, and the Isaac Newton 
Institute (EPSRC Grant EP/F005431/1) for support and hospitality via an 
EPDI fellowship, during which much of the present work was initiated. 

2 Identical particles in one and two dimensions 

We refer to the brief introduction in [38], the original references |54 [ 125 } [T9l 
155] . and the reviews [I3E71[MlllIlll3[231[12],fora complete introduction 
to the general concepts of identical particles and exchange statistics in one 
and two dimensions. Here we will jump directly to the consequences of the 
general theory (for scalar particles) outlined in the introduction, and make 
precise our mathematical assumptions and notation. 

As usual, we write H k (Q) for the Sobolev spaces of square- integrable 
functions on O C ~R n with square-integrable weak derivatives to order k. 
For a function space J- we generally write J~(X; Y) to emphasize that the 
functions map X — > Y, suppressing the latter argument whenever Y = C 
or if otherwise understood from the context. The space of smooth and 
compactly supported functions on Q is denoted by C^°(Q) and its closure 
inside H k (Q.) by Hq(Q,). Given a space J- of functions on the (traditional) 
iV-particle configuration space X = (K ) N , the subspace of functions which 
are completely symmetric resp. antisymmetric w.r.t. particle permutations 
will be denoted J-" sym resp. -Fasym- The domain of an operator A is denoted 
by®(A). 

2.1 One dimension 

Depending on which approach one takes to quantization [25j EE EH HI ED 
|2"71 12"2"1 HI] , identical particles in ID can be modeled as bosons, i.e. wave 
functions symmetric under the flip r i— > — r of any two relative particle 
coordinates r := Xj — x^, together with a local interaction potential, singular 
at r = and either of the form 5o(r) or 1/r 2 . We write (cp. e.g. [41]) 

V s (r):=2rj5 (r), V u {r) := ^ ~ 1} , (14) 

with statistics parameters r/,a G K, for the corresponding cases of Sch- 
rodinger- resp. Heisenberg-type quantization. These statistics potentials 
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should correspond to the choices of boundary conditions (in a sense to be 
made precise below) 

dtp 

— — = nib, at r = + , (15) 
or 

resp. 

ip(r)~r a , r^0+, (16) 

for the wave function ib at the boundary r = of the configuration space. 
Note that in this sense r/ = resp. a = represent bosons, while r/ = +oo 
resp. a = 1 represent fermions (in the bosonic representation, i.e. after 
factoring out the sign of the permutation). Suggested by such pairwise 
boundary conditions, one may define (cp. e.g. [H]) the total kinetic energy 
for a normalized wave function ib £ L s 2 ym := (g)J m L 2 (]R) of N identical 
particles on R to be Tg/ H := / k jv T^ K (ip;x) dx where 



T s %(V>;x):=I^-Vf + £ V s/H ( Xj - x k M\\ (17) 

j=l l<j<k<N 

with a corresponding kinetic energy operator 



1 N d 2 

j=l 3 l<j<k<N 

Formally, we need to specify domains i^s/H resp. iFg/H f° r the quadratic 
forms ip h-> Tgm(ip) and operators T s / H so that they are closed, respectively 
self-adjoint. Let (for arbitrary dimension <i) 

l\:={xe (R d ) N : 3 j / k s.t. atj = a; fe } (19) 

denote the diagonal set where any two particles meet, i.e. where the statistics 
potentials Vg/jj are singular. In the Schrodinger case, 

TsW = lf lWf^ + 27?/ |^| 2 rfx (20) 
2 Jrn 



i7 1 (]R Ar ) n L 2 ym , and hence also defines an associated self-adjoint operator 
Tg on L\ . Its domain is s.t. ib S @s is twice (weakly) differentiable on 



defines a closed and semibounded quadratic form on the domain f^g : 

-2 

''sym 

l N continuous at the diagonals A, and satisfies (cp. also [30] ) 
9 9 \ ( d d \ 

d^'dx-J ^=** ~{dx-~dx-J ^ =x * = mxj=Xk - (21) 

In the Heisenberg case, there is a complication in the choice of domains 
associated with the symmetry q4 1 — a. We therefore begin by considering 
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the quadratic form 




initially on the nice space of functions ij) G C^°(M. N \is) n T 2 ym . We can 
then use the substitution ijj = f a (j), /(x) := Y\j<k \ x j ~ x k\i an d a partial 
integration to write (see e.g. [37] ) 

/ |VV| 2 cfe= / fa(l-a)^P + a^) |^| 2 dx+ / |V<£| 2 / 2a dx, 
Jr^ Jr^ V J / / iR^ 

(23) 

where A/ = on R N \A and | V/|// 2 = 2 ^ i<fe |a?j - £fc|~ 2 - Hence we have 
the representation Tn(ip) = ^\\Qaip\\ 2 — for such if; and all a£R, where 
Q a denotes the vector- valued differential expression (cp. also [45 J) 

Qa ■= / a V/- Q = V-aW, Wfa) := J2( X J ~ a*)" 1 - 

Now, consider Q a : L 2 (R N ) -> V'(R N \£\;C N ) as a distribution-valued op- 
erator (note that since W ^ Lj oc (M ) we need to be careful and remove the 
diagonals), and define Th with domain := ^(Q™ ax ) to be the quadratic 
form Tn(ip) := ^WQ™*"^]] 2 associated to the maximal extension of the op- 
erator Q a on T 2 ym , i.e. 

Wr x ) := W G ^ 2 ym : QaV> G L 2 (R";C")}, 

with the identification L 2 (R W ; C N ) = L 2 (R N \i\;C N ). Note that this defini- 
tion of Th does not exclude the boundary behavior (fT6|) for a > —1/2, taking 
e.g. ip(x) = / a (x)e-W 2 G and that for a = 0, ^ H = H^R* \&) = 
H^ ym (R N ) (note that the symmetry requirement is important here). The 
self-adjoint operator Th := ^(Qo ax )*Qa ax is then defined with domain 

^H := W G ^ H : G ^((Q^T)}- 

Another option would be to define 2Th to be the Friedrichs extension, i.e. 
(Q™ m )*Q™ m , associated to the minimally extended operator Q™ in , its do- 
main ^(Q™ in ) being the closure of C™(R N \&) inside L 2 syra w.r.t. the 
form tp i — y WQaipW 2 . In this case we can explicitly characterize its do- 
main as @(QT n ) = ff 1 iSym (l Ar \A) for a / 1/2. This follows because 
for if) G C™{R N \£\) the identity ^3]) for a. + 1/2 and for a = 1/2 (the 
Hardy inequality) implies 

ClIIV^I < <C 2 ||W|| 
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for some constants c\,C2 > (depending ono/ 1/2), and hence taking the 
closure w.r.t. the i^-form is the same as w.r.t. the Q a -iorm. 

On the other hand, the following shows that the two extensions Q™ m 
and Q™ ax are the same for sufficiently large a (while they are in general not 
for small a, such as a = 0; see also 0]): 

Theorem 2. For a > 1 we have Q™ m = Q™ ax , with domain 

S) H = Hl sym (R N \&), 

and hence the operator Th = \ (Q™ ax )*(5™ ax is eguaZ to £/ie Friedrichs exten- 
sion T H = ^{Qa m TQa m - In the case N = 2 we have for all a > 1/2 i/taf 

gmin = Qmax ^ ^ = H^JR 2 \ A) . 

Proof. We claim that, for a > 1, i/) 6 must satisfy ^ E H 1 (R N ) as well 



as 

|2 



[ , — \ o < oo for all j ^ jfe, (24) 

y R iv (xj - x k y 

from which approximability in Hq(R n \£\)PiL 2 ym follows by taking ip e (x) := 
U j<k Mxj - x k )ip(x), where e (x) := 0(x/e) and € C°°(K;[0,1]) sym- 
metric and identically zero for \x\ < 1 and one for \x\ > 2. Then ip e E 
and 

- Q a i> = {\[4>e - i) + (vn^) ^, 

where ||(JJ <fi £ — l)ip\\ — > and \\(Y[ <f> e — ^)Qai>\\ — > as e — > by dominated 
convergence, and for the second term we have by §2£ 



<C - — — — dx -> 0. 

Hence, HQa^e — ^)|| — > as e — ^ 0. 

To prove the above claim, consider the pointwise a.e. in \£\ identity 

\Q a ip\ 2 = |W| 2 + a 2 W ■ W\ijj\ 2 - aW ■ ({Vip)ip + ipVip) 
= \Vip\ 2 + a 2 W 2 \ip\ 2 - aW • V\ip\ 2 

= \Vip\ 2 + a(a - l)W 2 \ip\ 2 - qV • (W\ip\ 2 ), (25) 

where in the second step we used that if) E H^ oc (R N \£S) since Wif> E 
L 2 oc (R N \A) and Q a ip = Vip - aWip E L 2 (R N ) (note that then 8. 



3 I 



(djV)V' - ^ WO = in D'(M \A) and V|V>| 2 E L 1 1 oc (R and in the 
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third step we have -V • W = W 2 = |V/|// 2 = 2J2 j<k \ x j ~ x k\~ 2 ■ Now, 
for any e > 0, let us define the wedge-shaped set 



T e := {x G R iV : x x + (N - l)s < x 2 + (N - 2)s < . . . < ijv-i + e < x N } 

N-l 

= P) {x G R N : Xj + e < x j+1 }, 

3=1 



with piecewise flat boundary dT £ = {jf = i{dT £ )jj + i 



.r- 



(dT e )jj + i := {x £ 1 : Xj + e = Xj+\ and Xk + e < a^fc+i for fc 7^ j}. 

By (|25p and the divergence theorem (see e.g. Theorem 6.9 in [31]) we have 
for e > 

|Q a V| 2 *c= / (\Vij\ 2 + a(a-l)W 2 \^\ 2 )dx-a [ h-W\4>\ 2 dx. 

(26) 

On the boundary component (dT e )j j+i, the unit outward normal is h = 
(ej — ej+i)/\/2, and we obtain pointwise 

v^n ■ W = - Wj+i 

iV 

= 2(xj - Xj+i)' 1 + ^ ((xj - Xk)' 1 - (x j+1 - Xk)' 1 ) 
k=l 

= -2s' 1 + ^2 t{xj - x k )~ l {xj - x k + a)' 1 
= —2e~ 1 +'^2 £ ( x j~ x k)~ 1 (xj—Xk+£)~ 1 + £ ( x k— x j)~ 1 ( x k— x j—^)~ 1 - 

k<j k>j+l 

For N = 2 we then have —an ■ W = a£ J 2, while for N > 2 the two 
additional positive sums above are bounded by 

j'-i i-i 



k=l k=l V ' 



respectively 

Ar N-j-1 



1 AT-j-l 



E ^-^r 1 ^-^-)- 1 -- 1 E k(kTT) =£ n-j 

k=j+2 k=l y 1 J 

and hence 

A Tir ae" 1 ( j-l N-j-l\ as' 1 N ae' 1 4 

-an-W > 2 - — _ , = > 



V2 V 3 N-j J V2j(N-j)~ V2 N : 
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which is strictly positive for all N > 2. Hence, since all terms in the r.h.s. 
of (|26p are nonnegative for a > 1, and the l.h.s. remains finite as e — > 0, we 
must have Vip, Wip E L 2 (To) (also in the case a = 1 we find Vip E L 2 (ro) 
and therefore E L 2 (ro) since QaV' G -^ 2 (ro))- By symmetry of ip, this 
also holds on the remaining parts of M N \£\ obtained by permutation of 
the coordinates, and this would then prove the claim if we only knew that 
Vip E T>'(WL ) is really a function on all of M. N and not also having some 
singular component supported on S\. The support must furthermore be of 
codimension one (cp. Lemma below). But this would require a codimen- 
sion one discontinuity of ip at one of the boundary components (STo) jj.fi, 
which is in contradiction with both Wip E L 2 (K N ) and the symmetry of ip. 
This proves the first part of the theorem. 

It remains to prove the stronger statement of the theorem for the special 
case N = 2 and a > 1/2 (this was also proved differently in [7J, in a one- 
particle form). We start by showing that \xi — X2\~^ip is (transversally) 
uniformly bounded if a > 1/2. Note that for x\ < X2, 

{x 2 - xi)" a V(x) = - / — ((s 2 - xi + 2t)~ a i;(x 1 -t,x 2 + t)) dt 
Jo dt 

/•oo 

= - (x 2 -x 1 +2t)- a ((d 2 - dx)ip( Xl -t,x 2 + t)- 2a(x 2 - x x + 2t)~ 1 ip) dt 
Jo 

/•oo 

= - (x 2 -x 1 + 2t)~ a ((Q a , 2 - Q a ,i)ip) (xi -t,x 2 + t) dt 
Jo 

1 

2 2 
\(Qa,2 ~ Qa,im (xi-t,X 2 +t)dt 



< 



(x 2 -xi + 2t)' 2a dt 







where we used Cauchy-Schwarz, and || • || r denotes the L 2 -norm w.r.t. the 
transversal variable r = X2 — x\ (but still depending on the longitudinal 
center-of-mass R = {X1+X2) /2). Similarly for x\ > X2, we have |r|~2-0(x) < 
||(Qa,2 - Qa,l)ip\\r/V^a - 2. Taking ip E (x) := (p E {xi - x 2 )^(x), with cp £ a 
logarithmic cut-off function defined as in the proof of Lemma [3] below, we 
have 

Qa.tpe ~ Qatp = {fe ~ tyQalp ± ^A- 

The norm-squared of the last term is 
.2 f f b'(elnr)| 2 



£ 




IrJt r 2 



drdR<Ce [ \\(Q a ,2-Qa,l)4>\\rdR 

Jr 

= Ce\\(Q a , 2 -Qa, 



which shows that \\Q a (ip e -tp)\\ -> as e -)• 0, with ip e E H% (M 2 □ 
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We will in the following always assume r/ > resp. a > 1 in order 
for the statistics potentials Vs resp. Vh to be nonnegative, and can hence 
work with the Friedrichs extension for Th- To emphasize the dependence 
on the statistics parameters, the corresponding domains will be denoted 
— -^s res P- — ^H- ^ * s well-known that in the exactly solvable 
Calogero-Sutherland models, involving the operator Tr, the eigenfunctions 
satisfy the boundary condition (fT6|) (see [3 [55]). However, for large a this 
requirement is too restrictive, since e.g. ip(x) = / 2 (x)e~' x ' is both in fPg 
and @S, even for a > 2. The explicit dependence on a of the domains 
and for the operators appearing in the two-particle case N = 2 have been 
thoroughly investigated in [7]. The complete behavior in the many-particle 
case for a £ (—1/2, 1) remains, to the best of our knowledge, an interesting 
open problem. 

2.2 Two dimensions 

(Abelian) anyons in R 2 comprise a continuous family of identical particles, 
characterized by a one-dimensional unitary representation of the braid group 
Bn, i-e. a complex phase e ia7T £ U(l), or a real statistics parameter «eR 
(modulo 2). For a = such particles are bosons while for a = 1 they are 
fermions. For general a, one can model quantum mechanical wave func- 
tions ip of N such anyons by means of bosonic wave functions on R 2 , i.e. 
completely symmetric functions in L 2 ((M. 2 ) N ), together with magnetic inter- 
action potentials Aj := a ^2k^j( x j ~ Xk)^ 1 ! of topological type between all 
particles. Here xl denotes the 90° counter-clockwise rotation of the vector 
x G R 2 , and x^ 1 := x/\x\ 2 . Hence, the total kinetic energy for ip is given 
by Ta := J R 2N T^(ijj; x) dx, where 

1 N 

T^(^x):=-J2\D^\ 2 , D j :=-iV j + A j 

3=1 

corresponding to a kinetic energy operator 

1 - 

T A-=^J2 D rDr (28) 

3=1 

There is also here a question on how this operator and quadratic form should 
be formally defined since Aj (x) diverges badly In [38j we defined 

the operator Ta in (|28p with a closed quadratic form Ta as the Friedrichs 
extension initially defined on C^°(M? N \i\)DL 2 ym . This might seem to imply 
a mild hard-core requirement for the wave function for o^O (cp. e.g. |36tl5] 
where a stronger form of hard-core requirement for anyons was imposed). 
Let us first observe explicitly that this is not the case whenever a £ 2Z 



iV Xj +a^2(xj-x k ) 1 I, 
(27) 



k=l 



14 



(respecting the periodicity of the statistics parameter), and then show for 
general a£l that this definition is truly non-restrictive in the sense that it 
coincides with that of the maximal extension for the quadratic form Ta. In 
other words, if we insist on the kinetic energy form (|27p being finite, then 
we necessarily have this self-adjoint extension. 

We begin by defining the magnetic derivative as the distribution-valued 
operator 

D : L 2 (R 2N ) V'(R 2N \A; C 2N ), jf, ^ (D^)f =1 , 

where Dj are the C 2 -vector- valued differential expressions in (|27p . Note 
that we again need to consider this as a distribution strictly outside the 
diagonals since Aj $ Lf oc {R 2N ) and hence Ajip <£ L} oc {R 2N ), although 
Aj e C°°(R 2N \A) and G V'(R 2N ). We have a corresponding max- 
imal domain for the operator D acting on the Hilbert space L 2 ym (R 2N ), 
which defines the form Ta: 

^A := ®{D™*) := ^ € L 2 ym : Di/> G L 2 (R 2N ; C 2N )}, T A ty) = 

This also defines an associated self-adjoint operator Ta := I(£) max )*£) max 
with domain 

#A := G S>k '■ Dip G ^((TJ max )*)}. 

To emphasize the dependence on a, the corresponding domains will be de- 
noted 

Again, another option is to consider the Priedrichs extension, i.e. define 
the minimal extension D mm as the closure of D acting on C£°(R 2N 

$(D min ) := C™(R 2N \A)nL s 2 ym C A , 

where the closure is taken w.r.t. the graph norm in L 2 yia x L 2 and the form 
V> ||T»V|| 2 - The associated self-adjoint operator is T^ 11 := ±(TJ min )*T> min . 
Note first that for a = 0, i.e. for the free kinetic energy for bosons, these 
two definitions coincide, Ta = T™ 111 = To, as a consequence of the following 
lemma. 

Lemma 3. We have equality for the Sobolev spaces ff 1 (R 2Ar ) = Hq (R 2N 
H^ ym (R 2N ) = Hl sym (R 2N \i,), Hi ym (R 2N ) = Hl asym (R 2N \is), as well as 
iJ 1 (R 2Ar ) = fl-^M^XA). 

Proo/. We first show that ^ Wasym) (M 2 ^) = ^ (sym/Mym) (M 2JV \A). Since 
C-(R 2Ar ) n L 2 sym/asym) is dense in n T 2 sym/asym) we can assume 

that ^ G C c °°(]R 2Ar ). Define ^ £ (x) := n j<fe y e (a!i - a>fcM*) for e > 0, where 
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<p e { x ) '■= \ x \)-> an d V 9 ^ C°°(]R; [0, 1]) is taken to be identically zero on 

(-00, -2) and one on (-1, 00). Then A G C™(R 2N \ A) D £ 2 sym/asym) , 

II& " V'll 2 = / in^ " 1| 2 IVf -> 0, e -> 0, 
since the measure of supp(n ip £ — 1) n supp ^ tends to zero, and 



l|v(v e - VOII = ll(n^ £ - i)vv + (vn^ £ 

< IKllVe - IJV^IUa^pp^) + ||Vn^ e ||L2(supp^)ll^l|L°°- 

The first term on the r.h.s. tends to zero as above, while for the second we 
have for each j G {1, . . . , N} 

\\VjUk<lM x k-Xi)\\ L2{supp ^ < \\Vj < Pe( x 3 -^lU^suppV) 



and 



iVjVe^-aJfcJIIia^upp^) < / / \V<p £ (x)\ 2 dxY[d Xl . 

J suppi/> J Mr 



— f \Vcp £ \ 2 dx = e 2 [°° \<p'(e\nr)\ 2 — = e 2 [ ^ \<p'(es)\ 2 ds < Ce, 

2tt Jr2 Jo r J-2/e 



'suppi/> . 

Finally note that, by the substitution r = e s 

-l/e 
-2/e 

which proves the first part of the lemma. 

It remains to show that we can identify H 1 ^™) = fl' 1 (R 2JV Any 
element of the l.h.s. can obviously also be interpreted uniquely as an ele- 
ment of the r.h.s. Conversely, assume that ip G L 2 (R 2N ) is s.t. the locally 
defined distribution Vip G T>'(R 2N \£±) is actually a function on R 2N \A 
which extends to be square- integrable on R 2N . Denote this function by 
H G L 2 (R 2N ). To see that this function is identical to the canonically defined 
distribution Vip G H~ 1 (R 2N ) when both are considered as distributions on 
the full space R 2N , let / := Vtp — fi G H~ 1 (R 2N ) and observe using the first 
part of the lemma that for any u G H 1 (®? N ) = H^(R 2N the dual space 
of fl- 1 ^ 2 ^), with u <- u n G C C °°(R 2JV 

(/,«) = lim (/,u n ) = 0, 

n— >oo 

since the distribution / is supported on A. Hence / = in T>'(R 2N ). □ 

Since for a = 0, by definition S>(D max ) = F^IR^ \ A) and S>{D min ) = 
H^ sym (R 2N we have by the above lemma that D min = D nmx and hence 
obtain the natural kinetic energy operator Tq = — for free bosons with 
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domain SPj^ = H 2 ym (R 2N ). Now, define the unitary multiplication opera- 
tor U : L 2 {R 2N ) -> L 2 (R 2N ), 

(C/V)(x) := ]J , Zj ~ Zk M x) = Yl e^^(x), x G M 2JV \A, 

. \Zj — Zk\ 

where := arg(^- — z k ), Zj := sc^i + ixj t 2- We then find that U : £ 2 ym — ^ 
-^asym) -^asym ~~ ^ -^sym i s an isomorphism and furthermore, noting that U 
is smooth on R 2N \A, we have d(<*=™+/9) = U- n D^ a =^U n as distribution- 
valued operators for n G Z and /3 G M. In particular, 

D (a=2n) = ^-2n(_ jv)z7 2n &nd f|a=2n) = {j-Znf^n^ 

with HD^ 2 ' 1 )^!! = ||V(£7 2Tl V)ll, and therefore 
Similarly, 

^-= 2n+1 = U^ 2n+1 ^Hl ym (R 2N ) and ®l= 2n+l = U^ 2n+1 ^ Hi ym (R 2N ). 

This also shows that we might have Vip $ L 2 (R 2N ) although Dtp G L 2 (R 2N ). 

The case of general a G R is a little trickier. Since Aj ^ L 2 oc (IR 2Ar ) we 
cannot directly apply the standard theorems for magnetic forms (see e.g. 
Theorems 7.21 and 7.22 in |31j). but they can nevertheless be extended to 
the anyonic case as follows. 

Lemma 4 (Diamagnetic inequality). Assume ip G L 2 (R 2N ) and Dip G 
L 2 (R 2Ar ;C 27V ) for arbitrary fixed a G M. Then \tp\ G H l {R 2N ), 

|V,H( X )|<|(^)(x)| (29) 
for all j G {1, . . . , N} and a.e. x G R 2N , and 

IIVl^HI^QfpiV.^jV) < \\Dtp\\ L2(R 2N. c 2Ny (30) 

Proof. Note that \ip\ G L 2 (R 2iV ) trivially. By the assumptions we have that 
Dtp, considered as a distribution acting on C^°(R 2N is actually a func- 
tion locally on R 2N \A and which is square-integrable on all of R 2N . In par- 
ticular, Dip = -iVip + Aipe L 2 oc (R 2N where also ,4^ 6 L 2 oc (R 2N 
and hence -iV^ G L 2 oc ( r27V \^)> Le - V> G ^(M 2 ^ By Theorem 6.17 
in [31] we then have \tp\ G fl'i c (M 2JV \A) and, for a.e. x G R 2N \A, 




Re(^w)(x), ifV(x)/0, 
0, if tp(x) = 0. 
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Now, for j G {1, . . . , N} and each such point x s.t. ^(x) ^ we have 
Re ( i- Vjt/i) = Re f 777 V.^ + iA^-fj-] = Re (ii(-%Vji/> + Ajip) ) , 



so 



-iVjip + Ajtpl, 



while for ^(x) = 0, (V|V>|)(x) = < \(-iVjil) + Ajip)(x)\. Hence, flU holds 
for a.e. x G R 2Ar \A (i.e. a.e. x G M 2Ar ), and 

|V i |-0|(x)| 2 (ix< f \(-iVjij + Ajij)(x)\ 2 dx. 

We have shown that G Hl oc (R 2N \A) with J r2 n \ V\ip\ \ 2 dx < 00, and 
hence \ip\ G ^(M 2 ^ \A) = ^(IR 2 ^) by LemmaEl □ 



Theorem 5. For arbitrary a £ R rae have that D mm = D max and hence 
f A = l(D max )*D max is equal to the Friedrichs extension f A = ±(D min )* D min . 

Proof. Take tp G Z$ A . We need to show that ip can be approximated by a 
sequence G C™{R 2N \A) nL 2 ym s.t. \\^ n -^\\ ^ and \\D{^ n -^)\\ ^0 
as n — )• 00. 

Step 1: We start by proving that we without loss of generality can assume 
ip G L 2 sym (lR 2Ar ), by finding a sequence ip n G L 2 sym (M 2Ar ) s.t. Dip n G 



L 2 (M 2Ar ), ||^ n - and ||I>(^„ - VOII ->• 0. We take := XnV' where 

Xjx) := x(x/n) and x G C c °°(M 2Ar ; [0, 1]) n L 2 ym with x = 1 for |x| < 1 and 
for |x| > 2. Then ijj n is symmetric with supp^n Q ^2n(0), 

HV'n-^ll 2 < / H 2 dx-^0, n^oo, 



as well as Dip n = -^Vx(x/n)V> + XuDip G L 2 (K ), and 

\W n - Dn < -II V X (x/n)^|| + ||(xn - 1 Wll < -ll^ll + \W\\V{B$ -> 0. 

Step 2: We next show that it is sufficient to assume ip G L^ sym (M. 2N ) . 
By step 1, we have i/) G L 2 sym (M 27V ) with G L 2 (M 2Ar ). Let ijj n (x) := 
XnOV 1 ! ( x ))V'(x) where Xn is the cut-off function defined analogously as above 
but on R. Then |^ n (x)| < 2n and ip n G L~(M 2Ar ). Furthermore, on R 2N \A 
we have ^ G i/f^ (M 2Ar \A) and, by Theorem 6.16 in [31] (cp. also Theorem 
7.22 in [E]), Vn G ^(M 2 ^ \A) and 

W„ = Xn(|V>l W + ^xM\)^\ E ^loc^ W 
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Hence, 



D^n=xM\)D^-i^M\)V\il>\ G L? oc (R 2iV 
eL°° eL°° 

with the full L 2 -norm of Dip n on M? N then bounded by 

\\XnM)\Um\\ 2 + HVx'ndV'DlloollVIVIIb < (1 + 2n)||£ty|| 2 < oo, 

where we have used the diamagnetic inequality, Lemma HJ Hence, ip n G 
and furthermore 

ll„/. „/.l|2 ^ I |„/.|2 



< / \ip\ 2 dx->0, 
J\ib\>n 



'\ip\>n 

and 

~ D^\\ < ||(X«(H) " + Ifo^M 

where the first term tends to zero by dominated convergence and for the sec- 
ond we have IV'IXnd^l) ^ C X{\tP\>ti}i hence also J R 2iv |^(I^I)| 2 |V|^|| dx ->■ 
0. 

Step 3: Our next step is to show that we can assume ip G L^ sym (R 2iV \£S), 
i.e. ip supported away from diagonals. By steps 1-2 we have ip G L^ sym (R 2iV ) 
and hence Dip G L 2 (R 2Ar ). As in the proof of Lemma [3] we define tpe{ x ) := 
Hj<k feixj - x k )ip(x) which is in L^ sym (R 2N \£S) for e > 0. We obtain 

= U<p e Di/> - t(V]> e )^ G L 2 C (R 2N 

and 

ll^n - < ll(n^ e - i)^ll + Il(vn^[|, 

where both terms tend to zero as e — > as in Lemma [31 

Step 4: By steps 1-3 we can assume ip G L^ sym (M 2Ar Hence t/> 
is compactly supported on some open set n := (R 27V \A) n B R (0), and 
also = -iVip + Aip G ^c(^)- We have then since G £ 2 (S1) that 
V> G i?o(^)> so there is a sequence V'n G C^°(0) s.t. 

W^n-iph^n) -*-0, and ||V^ n - VV>||z2(fi) -»■ 0, 

as n — )• oo. Hence £>Vn G L 2 (R 2Ar ) and 

H-Wn - V , )IIl2(m 2JV ) = l|-»V(^n - tp)+A(tp n - ^)|| L 2 (n) 

< ||V^„ - V^||i8(n) + IMU«(n)ll^n - V>IIl2(q) -> 0, 
which proves the theorem. □ 
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3 Local exclusion 



Since our wave functions are modeled using bosonic symmetry, but with 
the non-trivial exchange statistics represented by an interaction, we cannot 
take advantage of the usual exclusion principle encoded in the total (global) 
antisymmetry of the wave function. Instead we recall the following local 
consequence of the Pauli principle for fermions, given as Lemma 5 in |12j . 
and which was used by Dyson and Lenard in their proof of stability of matter 
in the bulk (see also [El [28]). 

Lemma 6 (Local exclusion for fermions in 3D). Let ip £ /\™L 2 (M 3 ) be a 
wave function of n fermions in M 3 and let Q be a ball of radius t. Then 

[ f^lV^I 2 ^ > (n-l)§ / |Vfdx, (31) 

where £ ~ 2.082 is the smallest positive root of the equation 

d 2 s'mx g 

dx 2 x 

We refer to such local energy bounds as a local exclusion principle since 
it implies that the energy is strictly positive whenever we have more than 
one particle, and hence that the particles cannot occupy the same single- 
particle state (which on a local region with free boundary conditions would 
be the zero-energy ground state). The inequality follows by expanding ip in 
the eigenfunctions of the Neumann Laplacian on Q, and it turns out to be 
sufficient with such a weak linear dependence on n. A corresponding family 
of bounds of this form for generalized exchange statistics forms the starting 
point for our proofs of kinetic energy and Lieb-Thirring inequalities for such 
statistics. 



3.1 Local exclusion for anyons and bounds for the energy of 
the ideal anyon gas 

The following local exclusion principle for anyons was proved in [38] , starting 
from a local pairwise relative magnetic Hardy inequality. The method then 
used for lifting such a pairwise energy bound to a bound for the full kinetic 
energy on a local region will be illustrated below for the case of intermediate 
statistics in ID (and is also in the case of fermions in three dimensions 
applicable for an alternative proof of Lemma [6] above, up to the value of the 
constant). 

Lemma 7 (Local exclusion for anyons). Let ip 6 be a wave function of 
n anyons and let SI C M 2 be either a disk or a square, with area Then 
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where cn is a constant which satisfies cn > 0.477 for the disk and > 0.358 
for the square, and C a ^ n is defined in ([6j). Hence, defining the local kinetic 
energy on Q for a normalized N -particle wave function tp E 

N 

T a:=E/ d^Vfxn(^)d Xj (34) 
we obtain the following local energy bound in terms of p 

^^wKi/^- 1 ),- (35) 

We shall consider an immediate application, partially along the lines of 
[12] , of this local exclusion principle for anyons in the form (|33p to an explicit 
lower bound for the ground state energy of a gas of non-interacting anyons. 
The numerical constant we obtain by this comparatively simple method is 
much better than the one following from the Lieb-Thirring inequality (|58p 
for anyons given below and proven in 



Theorem 8 (Ground state energy for N anyons in a box). Let ip E be a 
normalized N-anyon wave function supported on a square Q L CR 2 of side 
length L. Then 



i r N n 2 

T A = 2LEN 2& * C fc N Cl NJJ , (36) 



where cq is is the constant in Lemma [7] for the disk, and 

( 1 (1 + 2tA^\ 
cat := sup — 2 Tl, • 

7>0 7T 7 J 

Hence, the energy per unit area of the ideal anyon gas is bounded below by 

j2 > y CnC <*,nP ^ 0.059 C a N p, 

which, whenever the statistics parameter a = ^ is an odd numerator reduced 
fraction, remains finite in the thermodynamic limit and results in the bound 

rp -2 

-=£ > 0.059^, (37) 
L z u z 

as N — » oo and L — > oo while the total density p := N/L 2 is kept fixed. 
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Proof. Our approach is very similar to the proof of Theorem 8 in [12], which 
used the local exclusion principle (|3ip for fermions. Introduce an arbitrary 
length 1 > and write Q L j := Ql + B e (0). Then, 

N x 

2Ta = / 2N E l^l 2 ^2 / Xs/G/jte) dydx, 

j=l j Ql,i 

and after changing the order of integration and inserting the partition of 
unity 

N 

i = n (a:*) + XBKy) c = e n ^(y) (**) n xmvy o**) 

k=l AC{1,...,N} keA k^A 

(38) 

we find 



2Ta = ^L Y /^Ei^i 2 w*,) 

j Ql,i AQ {l,...,N} JR j=l 

X ]J XB t (y)( x k) ]J XB e (y)c( x k)dxdy 
keA k^A 

^2 j e / / e \°^\ 2 n ^ n 



'-^ AC{l,..,JV} (B<(y) : )JV _ |X | ^ keA H A 



We now apply f|33|) to each term in the first summation above, which involves 
a partition j4 of the N particles into n := \A\ of them being inside the disk 
Bi(y), while the remaining N — n residing outside, and therefore whose 
contributions to the magnetic potentials A.jeA can be gauged away. Thus, 
we find 



2T a>^I / E Cl^{\A\-l)_ 



AC{1,...,N} 

x J J IV'I 2 n dx k JJ dx k dy 



> / E l A l / 2N l^l 2 II XB t (y)(Xk) II XB e (yr(x k )dxdy 

j Ql,i AC{1,...,N} R keA k£A 

/ /,J^ 2 E II XB e (y)( x k) Y\ XB e {y) c ( x k) dxdy 



AC{1,...,N} keA k£A 
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We then revert the above procedure using (|38p and 

2 J^L H XB dy)^ X k)YlXB e {y)c{x k ) 
AC{l,...,N} keA k^A 

AT AT 

= Z ^(vJ^j) II Xs,(y)(a;fe) II Xb^vM 3 *) = 5^Xfl 4 (y)(a5j)» 

j=l AC{1,...,N} keA k£A j=l 



and hence 



2T A > 



> - + 2 ^) = cnCl N 9(^ 



setting £ := -^== with 7 > 0, and 



9(l) ■= —3. 



1 (1 + 27/VAO 



2 



7jvy2 7T 2 7 4 

This function has a positive maximum cat > 0.0023 for N > 2, while cat — >• 
1/4 as N —7- 00. In the case N = 1 we are left with the trivial bound on the 
energy. □ 

Note that for small N it might be preferable to use the Dirichlet bound 
for bosons instead (from the diamagnetic inequality (|29|) ) : 

V / |^| 2 dx>V/ \VM\ 2 dx>2N^[ |V| 2 ^x. 
, Jo? , Jo? L Jo? 



This resulting rough — but nevertheless non-trivial — bound (|37|) for 
the ideal anyon gas should perhaps be compared with the 'correct' ground 
state energy per unit volume for a non-interacting gas of density p of spinless 
fermions in two dimensions, which in our conventions is 

lim H = 7rp 2 . (39) 

Suggestions for improving the constant cq were given in [35]. Also note for 
comparison that the bound (|33[) for the disk holds with cq = 7r£ /2 ~ 10.65 
in the case that a = 1 (cp. with Lemma [6|), where £' denotes the first zero 
of the derivative of the Bessel function J\ . 
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Figure 1: Plot of Csiv) ( re d solid) and arctan y/rj + Arj 2 /ir 2 (green dashed) 
as a function of r] resp. t = — arctan(ry). 

3.2 Local exclusion in one dimension 

As for the anyons [38] we consider the local increase in energy due to the 
statistical interaction, starting from a pairwise relative consideration. 

Lemma 9 (Pairwise relative Schrodinger). Forr] > 0, consider the operator: 



d 2 

on a symmetric interval [—1, I] with Neumann boundary conditions atr = ±1. 
Then > ^sirf) 2 /I 2 > where £s(y) is defined as the smallest nonnegative 
solution to £tan£ = y given y > 0. 

Proof. Note that for rj > there are only non-negative eigenvalues A of H r] . 
The general solution to the corresponding boundary value problem 

-u" = Am, on (-Z,0)U (0,0, 
U '(Q+) = - u f (Q-) = r)u(0), 
u'(l) = u'(-l) = 0, 

with u(r) = u(—r) and A > 0, is u(r) = C ^cos y/Xr + sgn r sin \f\\ 

The condition at r = H demands ^/Xtany/Xl = rj, i.e. A = £,{i]l) 2 /l 2 is the 
smallest eigenvalue. □ 

Note that £s is monotone and £s(y) ~ y/y for small y, while £s(y) ~^ § 
as y — )• oo. A numerically very good approximation to £s for all y is (see 
Figure [Q) 

£ s (y) « arctan Jy + ^y 2 . (40) 
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We shall also find use for the following additional observations: 

Lemma 10. The map y i— > £s(y) 2 ; y > 0, is monotone and concave, so that 
for i] > 

%s{ri/x) 2 x > £s0?) 2 , a > 1, 

^{r]/x) 2 x < ts{ilf, 0<x<l. ' ' 

Furthermore, the map x i— > (,s(v / x ) 2 x3 j x >0, is monotone and convex for 
any rj > 0. 

Proof. Differentiating the identity £tan£ = y w.r.t. y > we obtain £'(£ 2 + 
y(y + 1)) = £ 3 hence £'(y) > 0- Differentiating once more gives 

e ,/ (e 2 + y(y + i)) = -2e , (^' + y), 

so 

(e 2 )" = 2^ + 2££" = ; ^ ; - f -= ; ^ ; - + y - ^ < 0, 



e 2 + y(y + i) \e + y(y + i) 2 

since the expression in brackets is positive for y > 0. The estimates (|4ip 
then follow since £(0) = 0. The proof of convexity for £(r//x) 2 x 3 follows by 
similar analysis. □ 

Lemma 11 (Pairwise relative Heisenberg). For a > 0, consider nonnegative 
realizations of the operator: 

d 2 a(a - 1) 



dr 1 r 2 

on a symmetric interval [—1,1] with Neumann boundary conditions at r = 
±1. (In the case —1/2 < a < 3/2 one has to pick the correct self-adjoint 
extension, matching the b.c. ()16p — see below, and Section \2.1\ ) Then 
H a > Cn( a ) 2 A 2 ; where £h(o0 is defined as the smallest nonnegative solution 
to 

J(0 + 2£J'(£) = 0, 
where J = J a _i is the Bessel function of order a — \. 

Proof. The general solution to the corresponding boundary value problem 
on [0,1]: 

„ a(a — 1) 
— u H g u = on (0, Z) , 

u'(i) = 0, 

is via the ansatz u := r?v found to be (we have assumed A > and focus 
on the non-degenerate case a / 1/2) 

u(r) = Cir5 J a _i (y/\r) + C 2 r^ Ji_ a (V\r). 
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Figure 2: Plot of £h(c*0 (red solid) and y / 7t 2 /4 + a (a — 1) (green dashed) 
as a function of a > 0. 

For a > 3/2, 6*2 = is enforced by the requirement of square-integrability, 
while for a < 3/2 this choice is taken to define the domain of H a . Note that 
for a > 1/2 the choice C2 = is required in order for r a (r~ a u)' £ L 2 ([0, 1]) 
and hence for (a corresponding two-particle version of) u to be in the domain 
of Q™ ax . The Neumann condition at r = I then demands J a _i(V\l) + 

2y/XU' _ L (y/Xl) = 0, i.e. A = £(a) 2 // 2 is the smallest eigenvalue. □ 



Note that according to the above £h(0) = 0, while £h(1) = f ■ A numer- 
ically very good approximation to £h on the interval [0, 1] is 



However, note that £h continues to grow approximately like £h(«) ~ a for 
a > 1 (see Figure [2]) . In any case we have from the lowest eigenvalue of the 
Laplacian with Dirichlet boundary condition at r = 0, Neumann at r = I, 
and with a potential > a(a — l)/l 2 '- 



Now, consider the Neumann kinetic energy for n > 2 particles on an 




< a < 1. 



(42) 




for a > 1. 



(43) 
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interval Q = [a, b\. Using the identity 



nY,\dM 2 = E \d^-d k ^\ 2 + 

j=l ^<j<k<n 



E^ 

i=i 



(44) 



to separate out the center-of-mass motion, we find 



/ 

JQ 



Tg /K (ij;x)dx> 



[ ( ^ E - d ^ 2 + E f s/h(^- - x k m 2 ] 

QU \ j<k j<k J 



f/x 



n ^ 



-' ./Q J[-S(R),S(R)] 
2 / / <-S/H 



(|a r vl 2 + u s/ H(OIV'l 2 )dr ( ii^ ^ ix , 



^-E I I 3^2 I M 2 drdRdx>, (45) 

n f^ k JQ"- 2 JQ 6 ( R ^ A-S(R),5(R)] 



where we for each pair (j, k) split the parameterization of Q n into n — 2 
variables x' = (xi, . . . . . . . . . , xjv), a pairwise center-of-mass coor- 
dinate R := (a^ + £fc)/2 € (J, and a relative coordinate r := (cc.,- — x k ) € 
<5(i2)]. We have set 5(R) := 2min{|i?-o|, \R-b\}, and in the above 
we also used that Vg m > followed by Lemma [9] or QTJ 

In the Heisenberg case £h = £h(c>0 is constant on the interval Q and a 
local exclusion principle follows in a similar way as for the anyons: 



Lemma 12 (Local exclusion for Heisenberg). We have for a > 1 

,£ H (a) 2 



/ T£(^;x)dx > (n-l)^g- / 



(46) 



and hence, defining the local kinetic energy on Q for a normalized N -particle 
wave function ifi 6 



N 



\ d M 2 + E V H(Xj-%k 

0V)fc=i 



XQ(x j )(ix, (47) 



we obtain the local energy bound 

£h(«) 2 



rpQ \ 

J H - 



\Q\'< 



p{x)dx 



1 



Q 



(48) 



Proof. We use 5(R)~ 2 > \Q\~ 2 in (gSJ) to obtain (|4"6"j) . Then we proceed 
as for Lemma 8 in [38] (cp. also |12] and the proof of Theorem [8] above), 
inserting the partition of unity 



AC{1,...,N} l&A l$A 



(49) 
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in (g7D: 



A jR j&A \ U¥=)k=i j leA if A 



e H (g) 

IQI 2 



t(E^)-i) 



2 (fx, 



hence (|I8|) . where we again used Vh > and in the last step the partition 
of unity (gH). □ 



In the Schrodinger case, we observe from LemmaflOlwith x := \Q\/5{R) > 
1 that 

UvS(R)) 2 . ts(v\Q\/x) 2 x t s ( V \Q\y 

m 2 - — w — ~^qF~ foralli?G ^ ^ 

and can thereafter proceed exactly as above for the Heisenberg case. 
Lemma 13 (Local exclusion for Schrodinger). We have for r\ > 



T s "(^;x)dx > (n-l) gs( ^ l)2 / jrj-r/x. (51.) 



and hence, defining the local kinetic energy on Q for a normalized N -particle 
wave function tp E f^<? 



3=1 M \ (?¥)*=! 



(a?j -xjfe)|^| |xQ(zj)dx, (52) 

we obtain the local energy bound 



Zs(ji\Q\) 2 



T g > «*^H / -1) . (53) 



Proof. Use the estimate (|50p in (|45f) and continue as in the proof of Lemma 
E2J □ 
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4 Local uncertainty 

We shall use the following local form of the uncertainty principle on a d- 
dimensional cube Q. 

Lemma 14 (Local uncertainty principle). For any e G (0, 1) we have 

T Q > 9k £ ^ - S (i + ( -EJ) 1+4/ ^ V , (54) 



where the local kinetic energy denotes either with rj > 0, or lozi/i 
a > 1, or wt/i any a G R, or simply Tq , the free kinetic energy on a 
d-cube Q in any dimension d > 1 (i.e. half of the l.h.s. of (I55p ). 

This follows from the following lemma concerning the free kinetic energy 
for an arbitrary iV-particle wave function, which was proved as Theorem 14 
in [38] using methods from |48| . 

Lemma 15 (Uncertainty on a cube). Let Q be a cube in M. d with volume 
\Q\, and let u G H 1 ^^) be an arbitrary N -particle wave function. Then 

i 2(d+2) 

eLiv^ 'far 



p(xy- 



\Q\ 



d 

dx, 



(55) 

where p(x) := J2j=i Ir^n-i) \u(xi, . . . , Xj_i, x, x j+1 , . . . , £Cat)| 2 Y[ k ft dx k- 

_2 

The constant C' d := d 2 C d d / (d + 2)(d + 4), with Cd given below, enters 
via a bound for the growth of the sequence of eigenvalues {Afc}£l for the 
Neumann Laplacian on Q: 

?L<C d ei, e>0. (56) 

0<A fe <e M 

In the case d = 1 we have A^ = j^k 2 , k = 0, 1, 2, . . . , and 

#{k : < A fc < e} = #{k G Z : < k < e2|Q|/vr} < e^\Q\/ir, 

and hence we take C\ := 2/tt and C[ = 7r 2 /60. In the case d > 2 we have 
Afc = 7r 2 |fc| 2 /|(5| 3 , k G Z> , and by geometric considerations 

#{fc : < A fc < e} = G Z> : < |fe| < e3|Q|aM < d(e* |Q| s/tt) ', 
i.e. the optimal constant in (fH>6"j) is := d2 d /ir d , so 

^ = T(d + 2)(d + 4)' (5?) 

In particular, C2 = 8/-7T 2 and C 2 = it 2 / 48. We remark that these constants 
are probably far from the sharp ones in the inequality (|55p . 
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Proof of Lemma [14\ We reduce to the case of free bosons, T Q > T Q , by dis- 
carding the statistics potentials or using the diamagnetic inequality (Lemma 
HI \Djip\ > (VjlV'll for anyons (which also preserves p). Then, from Lemma 

m 



J o - 



2(foP) 2/d 



p(x 



1-1 2+4/d 
2 



\Q\ 



dx, 



with 



2+4/d 




1+2/d 



dx 



UoP) 



l+2/d 



2+4/d 



\Q\ 2 / d 



where the norm 



i 

\P 2 \ 



2+4/d 



\\(I Q p/\Q\) 



p is that of L P (Q). The first term is bounded below by 

2+4/d 



\2+4/d 



using the triangle inequality. Fur- 
e G (0, 1), and p > 1 



thermore, by convexity we have for any a, b G 
that 

(ea + (1 - e)b) p < eaP + (1 - e)P, 
and hence with a = A — B and b = jt^B, 

(A - B) p > e v ^ x A v - (^Y~~^j P BP - 

Applying this inequality to the norms above with p = 2 + 4/d, we finally 
arrive at (15H). □ 



We remark that bounds of the form (|54p could also be obtained by 
application of standard Sobolev and Poincare inequalities, but would then 
need to be dealt differently for d = 1,2 and d > 3 (see [2] concerning the 
latter case). For d = 2 an alternative form of the local uncertainty principle 
than (j54"|) was given and used in 



5 Lieb-Thirring inequalities 

5.1 A Lieb-Thirring inequality for anyons 

In [38] we introduced an approach to Lieb-Thirring inequalities based on 
local exclusion and uncertainty, resulting in the following energy bounds for 
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anyons in M 2 . 



Theorem 16 (Kinetic energy inequality for anyons). Let ip £ St^, a € M, 
be a normalized N-anyon wave function on M. 2 . Then 

[ if (^;x)cfcc > C A C 2 aN [ p(x) 2 dx, (58) 
Jr 2N Jr 2 

for some positive constant 10~ 4 < Ca < vr. 

Corollary (Lieb-Thirring inequality for anyons). Let if) £ S> A be a normal- 
ized N-anyon wave function and V a real-valued potential on M 2 . Then 

[ \Tg(t,x)+J2 V ( x iM 2 | > -C' A C~ 2 N [ \V-(x)\ 2 da, (59) 
Jm. 2N \ ~[ J Jr 2 



for a positive constant C' A = (4Ca) 



-i 



The rough lower bound Ca > 10 -4 that is given here for the optimal 
constant Ca in f)58[) can be computed in a similar way as we do below for 
the Schrodinger case (in the proof of the corresponding Theorem 11 in 
one can take e.g. e = 1/8 and k = 7cqC 2 n /(Att 2 )). Note that from ([5? 
follows also the bound for the free anyon gas (see [38] ) 

p > C A Cl N p 2 , (60) 

with mean density p = N/\Q\, assuming that p is supported on a domain 
ft C M 2 . Comparing with (|39[) it is clear that the optimal constant in (j58[) 
and (|60p satisfies Ca < C A := ir, which is the exact constant for fermions in 
the semi-classical approximation. We should also compare with Theorem [5] 
which produced a significantly better numerical bound for the constant for 
the anyon gas. 

5.2 Lieb-Thirring inequalities in one dimension 

By proceeding as for the anyons, combining local uncertainty with local 
exclusion, we can prove the following for intermediate statistics in one di- 
mension, summarized in Theorem [TJ 

Theorem 17 (Lieb-Thirring inequalities for ID Schrodinger). For tp G 3)^ 
with rj > and N > 1 



7f ty;x)cfc > C s / ^(2nrp(x)?p{xfdx, (61) 

JV 



for some positive constant 3-10 5 < Cs < 2/3, where p is a piecewise 
constant approximation of p, defined below via a splitting algorithm. This p 
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in ()6ip can also be replaced by the Hardy-Littlewood maximal function of p. 
In particular, if -rip- < jp for all intervals Q and some 7 > 0, then 



I 

Jr 



2^;x)cfcc > C s es(2r?/(7P)) 2 / p(xf dx, (62) 

Jr 



Jr 1 



N 



(63) 

with C g := |(3Cs) _ 2, and i/ p is supported on an interval of length L 

T s /L > C s ^s(2r ? /( 7 p)) 2 /5 3 , p:=N/L. (64) 

Proof. For N = 1 and iV = 2 we use Ts > To and Lemma [15] (which 
generalizes to Q = K) with a resulting constant C{/8 > 0.020. 

For N > 3 we can consider the kinetic energy on an arbitrary finite 
interval Qo C K s.t. Jq q p > 2. We split the interval Qq in halves iteratively, 
organizing the resulting subintervals Q in a (full binary) tree T; cp. |38| . 
The procedure can be arranged so that Qq is finally covered by intervals Q b 
marked B s.t. 2 < Jq p < 4, and marked A s.t. < Jq p < 2, sitting 
at the leaves of the tree and s.t. at least one B-interval is at the highest 
level of every branch of the tree. On the B-intervals we use local exclusion, 
(|53|) . together with local uncertainty, Lemma [21 to obtain 

C[Jq b P 3 C' x (,( e \ 5 \ 4 , 1 M(v\Qb\? 



(65) 

for any e, k G (0, 1). For simplicity we set e = 5 and « = ^s( ? ?|QB|) 2 /( 7r /2) 2 
< i to find 



■ c;/2 10 = 2- 11 / 

Note by monotonicity of £s that 



with ci = 2/vr 2 • CJ/2 10 = 2- n /15 and c 2 = -2/vr 2 • 4C[ + 1/2 = 11/30. 



r 

Cs(r/|Q B |)>es(2n/p(x)), where p| Qfl := 

i.e. /5 is defined to be the mean of p on each B-interval. The A-intervals are 
further divided into a subclass A2 on which 

p 6 > c- * 



\Qa\ 
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with c := 2 7 , so that by Lemma [Lil (again using e = 5) 

Tq A2 > c 3 / P 3 = ^ [ Zs&v/p)V, where j5| := 0, (67) 
JQa 2 77 JQa 2 

C3 = ^C[/2 8 = 7r 2 ci, and a remaining subclass Ai for which 



/ 

JO 



A 



„3 ^ „ vj Qa 

|Q 



Consider the set A\{Qb) of such intervals Qa 1 which can be found by going 
back in the tree T from a fixed B-interval Qb at level k £~N (possibly all the 
way to Qo), and then one step forward. On each level 1 < j < k there is at 
most one such interval which we denote by Qj (and otherwise we can define 
\Qj\ := IQol/2- 7 for the below expressions to make sense). Using Lemma [TPl 
we have 



2 \Qb\ c , lr , ,x2 



e s (ji\QB\/(\QB\/\Qj\)y ^ 



and hence 



j=l Ml j=l \^3\ j=1 

< 2Z S ( V \Q B \) 2 . 

Defining p\q Ai ■= Jq a p/\Qai I to be the mean also on Ai-intervals, we have 
by < p < 2/IQxil an d monotonicity of £s(rj / 'x) 2 'x 3 that 

/ &(2 V /p) 2 p 3 < ^(2 V /p) 2 cp 3 \Q Al \ < c^^lQ^I) 2 -^-, 
JQ Al MM I 

and therefore by the above 

for every B-interval Qb- In other words the energy on all intervals with 
almost constant low density is dominated by that from exclusion on the 
B-intervals. Using (j66|) we then have 



T^> Cl [ Cs(2r ? /p)V + c 2 £ / es(2r//p) 2 p 3 , (68) 



c ' 2 = c 2 /(16c) = 2~ 12 • 11/15, and hence by flBT]) and ([68]) 

T S > T s Qo = £ T^ B > C S f H2n/~P?P\ 
Q A „eT J Q° 



33 



where Cs '■= min{ci, c' 2 , Acs/tt 2 } = c\ > 3- 10~ 5 . Now take the limit \Qo\ — > 
oo and arrive at ()6ip . 

The bounds (f63j) and (f64"P follow in a standard way from (f62j) (see e.g. 
Theorem 4.3 in [33]). The bound Cs < 2/3 for the optimal constant in 
(|6ip follows by taking the limit r/ —> +oo and comparing (|64p with the 
semiclassics for fermions. □ 

Theorem 18 (Lieb-Thirring inequalities for ID Heisenberg). Let ijj £ f^-g 
mt/i a > 1 and N > 2. Given any finite interval Qo s.t. Jq q p>2, we have 

Tg° > C H £ H (a) 2 / ~p{x?dx > CMaf {!Qo ^ X) \ 
jQo \Qo\ 



(69) 

for some positive constant 1/32 < Ch < 2/3 ; and p defined below (see (|72p ) 
as a piecewise constant approximation to p. In particular, 

Tn = I Tg(i;;x)dx > C H £ H (a) 2 [ p(xf dx, (70) 

JR N JR 

and if p is confined to a length L 

T H /L > C H £ H (a) 2 p 3 , p:=N/L. (71) 

Here the assumption N > 1 is important because of the dependence on 
a in the bound (|70|) . Note that the integral in the r.h.s. of (|70|) actually 
becomes a Riemann sum approximation of J* K p(x)^ dx in the limit of many 
particles with relatively uniform density. Also note that the usual Lieb- 
Thirring inequality ([2]) for fermions (with a constant independent of a) is 
valid as a lower bound for Th for all N and a > 1. However, we will prove 
below that p cannot be replaced with the actual density p in the r.h.s. of 
([70j) as long as Ch is a positive constant independent of a, using that £h(o) 
can become arbitrarily large with a — > oo. 

Proof We split an arbitrary finite interval Qo s.t. p > 2 iteratively as in 
the proof of the previous theorem (in the case N = 2we can approximate this 
condition arbitrarily well with a finite interval, and some minor adjustments 
to the constants below will need to be incorporated accordingly). In the 
Heisenberg case we cannot combine (|48|) with Lemma [TH to yield a uniform 
bound w.r.t. a (better than the one for fermions) since £h(c>0 — s- oo as 
a — > oo. Instead we need to rely solely on energy from exclusion (also 
resulting in a slightly better lower bound for the constant Ch). We define 

p1q TqT (72) 
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for all A- and B-intervals Q G T, and use for < k < 1 

(73) 

(with x — 1 > 3x 3 /4 3 when 2 < x := Jq p < 4) for the B-intervals, and 



'' (/<?i^ 3 8 8 4 fc+1 2 5 /3 



< 



2 



for all the A-intervals associated to a B-interval Qb at level in T (similarly 
as for A\{Qb) in the previous proof). Hence, 



2 5 — 




and 

T H > > ^ > C H £ H (a) 2 / p 3 , 

Q B eT 



where « := 2/3, C H := 2~ 5 . 

For ([MI) and flUi we use that f Q p 3 > \Qo\~ 2 (Jq p) 3 and f Q p = f Q p. 
The upper bound on the optimal constant Ch follows by comparing (|7ip 
with a = 1 to the semiclassical limit for fermions. □ 

Theorem 19 (Necessity of p in (|70p ). Assume that for some nonnegative 
constant Ch f/ie inequality 



[ r^(^;x)dx > C H £ H (a) 2 / P(xfdx 

JR N JR 



(74) 



holds for all ip £ and a > 1. TTien we must /ia?;e Ch = 0. 

Proof. The idea is to spread out the particles one by one but with high 
individual localization, and then compare the behavior of the left and right 
hand sides of (|74p as a — > oo. Let (p £ C£°([— 1, 1]; R>o) be a bump function 
normalized s.t. J^_ l ip 2 dx = 1, and define cp E (x) := e~2(p(x/e), and the 
completely symmetric iV-particle wave function 

1 n 

ViV - <res Jv i=i 

Assuming < £ < 1/3 we have on the support of tp that \xj — Xk\ > 1/3 
for all j k, so this wave function is certainly in the domain of Th- It is 
furthermore normalized: 

1 /■ n 

l^l 2 dx= l^ Yl / IT Ve{ x 3 ~ v(j))Ve(Xj ~ T(j)) dx = 1, 

Nl „ Z7i Ju N 7*i* » ' 
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with one-particle density 



p(x) = N I \ip(x,x 2 ,...,x N )\ 2 dx' 

jJV-l 



N f N 

N - v,T£S N J^- 1 3=2 

N 

= ^2 ( Pe(x- jf, 

i=i 

supported on the interval [Q,N + 1], and s.t. p(x) 3 = YljLiVei 31 — i) 6 - 
Furthermore, 



and ^ 

= S e ~ V((a:j - II Vs(xj - <r(J)), 

implying J rN \dkip\ 2 dx = e~ 2 f\(<p') 2 dx. Hence, 

T u <e- 2 N [ {ip') 2 dx + -a(a- l)iV(iV-l), (75) 
J-l 2 



while the r.h.s. of (|74p is 

JV 



C H ?h(«) 2 V / e~\{{x - j)/ef dx = £- 2 NC n £ H (a) 2 / if 6 dx. 
~[ J-i 



(76) 



Now, taking e depending on N so small that e~ 2 NC K (p e dx > 9N(N-1), 
we find using (|4"3"j) . £h(") 2 > a (a — 1), that (f76|) is strictly greater than (f75|) 
for a sufficiently large, unless Ch = 0. 

Note on the other hand that p, and hence the r.h.s. of (|70p . is indepen- 
dent of e for the above choice of ^. □ 



6 Some applications 

6.1 Many anyons in a harmonic oscillator potential 

As an application of the kinetic energy inequality (|58|) for anyons, we can 
consider the ground state energy of N anyons in a harmonic oscillator po- 
tential, a problem which has been discussed substantially in the literature. 
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The Hamiltonian operator is in our conventions given by 



N /I \ 

.7=1 V J 



V(x) :=^\x\ 2 . 



For N = 2 the exact spectrum was computed already in [25], while for iV > 3 
the ground state energy is only known exactly for fermions a = 1 and in a 
small neighborhood around bosons a = 0. The lower part of the spectrum 
has been computed numerically for the full range of a G [0, 1] for N = 3, 4 
[52|, |4U1 [53] and for general TV" there are some known families of exact but 
excited or singular eigenstates [5"8 | [TD | 15]. Some general features extend to 
arbitrary N, such as large numbers of level crossings in the ground state as 
TV" is fixed large and a varied [9] . See [23] for a more recent summary of the 
status of this problem. 

With the kinetic energy inequality f)58[) at hand, we would like to find a 
lower bound for the quadratic form 

> j^(c A Cl N p(x) 2 + ^-\x\ 2 p(x)]dx (77) 
subject to the conditions J K 2 p = N and p > 0. We extremize the functional 



F[p,X] :-- 



J^ d (c A Cl N p{x) 2 + ^-\x\ 2 p(x) - Xp{x)\ dx + XN 



to find 



2C A Cl N p(x) + ^-\x\ 2 - A = 0, 



and 



N 



Jr 2 



x) dx, 



i.e. 



with 



2C A C 2 aiN N = 2vr 



p(x) 



[X-u 2 \x\ 2 /2\. 
2CaC 2 n 



2\/uj 



(A - u r /2)rdr = 2vr 



2\/w 



ttA 2 



and hence A = u)C a nW2C a N/tt. Plugging in this form for p, the r.h.s. of 
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dZZD is 



2tt 



2X/uj /i 



AT JO 
7T 



2 (A - oj 2 r 2 /2) 2 + ^-r 2 (A - u; V /2) ) r 



2 2 + 4 



2 



2A/w 



7T 



2C A C 2 



iV 



A^ 



3w 2 



7T 



(u;C Q , A r V / 2C A A70F) J 



3C A G 



a,7V 



and hence 



C~ N uN*. 



7T 



(78 



This resulting bound for the energy in terms of C a ^ujN^ allows us to com- 
pare the graph of a \-t C aj N (the limiting graph for C a := limAr_ s>00 C a< N is 
sketched in Figure 2 in [38]) with the previously found exact and numerical 
spectra for 2-4 anyons and the extensions of certain features in such spectra 
to more anyons (note e.g. the partial reflection symmetry about a = 1/2; 
cp. [12]). The correct ground state energy in the case a = 1 is to leading 
order ~ ^-ujN2. We also note that for bosons, which satisfy the weaker 
kinetic energy inequality Tq > J r2 p 2 , the above bound reduces to 



The correct ground state energy for bosons is of course ojN . 

The bound (|78p improves, for odd numerator rational a, the previously 
known best lower bound for the energy of ./V anyons in a harmonic oscillator 
@9l|9]: 

N(N - 1) 



{il>,Hi/)) > u [N + 



L + a- 



(79) 



where L is the total angular momentum of the state ip; (— i ^ • Xj A Vj)?/> = 
Lip. We see from this inequality that the energy actually grows like N 2 , 
unless L ~ — a( 2 ) (equality is possible for certain N and a) for which ([79]) 
reduces to the bosonic bound for the energy. The non-trivial dependence 
on a and N of these energy bounds, together with other aspects of the 
spectra mentioned above, raises the question whether the limiting constant 
in the kinetic energy inequality (|58p actually cannot be improved for even 
numerator and irrational values of a (for which C Qj at — > as N — > oo). We 
refer to [39] for further discussion on this possibility. 
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6.2 Stability of odd-fractional anyonic matter with Coulomb 
interactions 



In the context of the fractional quantum Hall effect it is relevant to con- 
sider anyons arising in the form of quasiparticles with fractional charge and 
statistics [2] embedded in a three-dimensional system, possibly with an effec- 
tive 3D Coulomb repulsion. Since these are then arising in a material there 
could also be a large number of oppositely charged particles which impose a 
Coulomb attraction on the anyons. We can ask if such an interacting system 
is stable in the thermodynamic limit, i.e. if the negative ground state energy 
grows at most linearly with the total number of particles (cp. stability for 
ordinary fermionic matter [33]). Consider therefore the model Hamiltonian 



1 - 

H(N, R) := — V Dj ■ Dj + V c (x, R), 
2m ^-^ 



z ^ z 2 



l<j<k<N 1 J 1 j=l k=l 1 J K| l<j<k<K 1 3 K| 

which can be taken to describe N anyons with mass m and unit charge — 1 
moving in the plane M 2 , together with K static 'nuclei' of charge Z > 1 at 
positions R = (J2i, . . . , Rk) £ M 2X , and where all particles are interacting 
through electromagnetic forces in 3D. The condition Z > 1 is only technical 
and can be relaxed, while the overall scale of the unit charge can be adjusted 
with the variable mass m. 

By the usual electromagnetic screening in 3D, we have the following 
important lemma to simplify the problem (given e.g. as Theorem 5.4 in 



Lemma 20 (Baxter's electrostatic inequality). 

where := ^ mm,yfc \Rj — Rk\, and dn(x) := minfc = i i ... j x \x — Rk\. 

Hence, we would like to find a lower bound, independent of R, for the 
quadratic form 

>—[ Y(\D J i;\ 2 -2m'^^\i;\ 2 ) dx+—Y-±-. (80) 
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With the Lieb-Thirring inequality (|59p for anyons at hand, we will prove 
the following theorem. 

Theorem 21 (Stability of 'anyonic matter' for odd- fractional statistics). 
For a normalized N- any on wave function -ip £ with a = fi/v an odd 

numerator reduced fraction, we have the lower bound 

(ip,H(N, R)V>) > -Cv 2 mZ 2 {K + N), 

for some positive constant C. 

There are some technical difficulties related to the dimensionality. In 
order to handle the non-square-integrability of the one-particle potential in 
(|80p close to the nuclei, we split the integral into two regions. On a region 
Be '■= UfcLi (-Rfc)i < £k < Dk, close to the nuclei we use \Djip\ > 
|Vj|^|| together with the following lemma to deduce 

f ( 1 1Z + 1 \ 

/ -\D 3 i,\ 2 - 2m—-^\^\ 2 XbAvj) 

>-2m(2Z + l)f (Am(2Z + l) + ^-]^\ 2 XBAxj)dx, (81) 

where e(x) := et on the disk B ek {Rk) around the kth. nucleus. 
Lemma 22. For v G H 1 (B £ (0)) and e,/x > we have 

I ( ^\Vv\ 2 - -^-\v\ 2 ) dr > - ( - + - J / H 2 (fr. 

Js £ (o) V M / \/* e/ ./b,(o) 

Proof. We simply adapt the proof of Lemma 2 in [12] to two dimensions. It 
is clear that the ratio 

( fi\Vv\ 2 - -^-\v\ 2 ) dr / [ \v\ 2 dr 

B s (0) \ \ r \ ) I JB E {0) 

is minimized by a radially symmetric function v(r), and that we may con- 
sider a lower bound when the domain is replaced by the annulus B £ (0)\B$(Q) 
if the bound is uniform in 5 > 0. We have therefore reduced to the Neumann 
problem 



-fi[v" + -)-- = Xv, v'(S) = v'(s) = 0, 
\ r J r 

on [5, e] , and would like to show that the lowest eigenvalue A > — - — | . 
Noting that the ground state v is nodeless and defining u) := —v'/v, we have 

-d r (ru) - fiu 2 - - = A, u(S) = u(e) = 0, (82) 
r r 
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and after multiplying by r and integrating over [5, e] 



— fj, / uj 2 rdr — (e — 5) = X- 



: 2 -5 2 



(83) 



Now, differentiating the differential equation in (J82J) by r we have 

■ 0. 



1 \ 1/1 

2u + 



r 2 \ fj, 



and therefore an extremal point ro £ (<5, e) for w cannot be a minimum if 
(jj{tq) < and cannot be a maximum if w(ro) > 1//U, hence < w < l//i. It 
then follows by (j83|) that A > — •= — which proves the lemma. □ 

Now, combining the integral term in (|80p with (|8ip using half of the 
kinetic energy, produces the lower bound 

(ip,H(N, R)ip) > 

2^ ^ 2JV E Ql^l 2 " 2m ( 2Z + X ) (2(2m(2Z + 1) + e-^XB.fo) 



> 



L f y ( -\DM 2 - V x (x i 



with the one-particle potential V\(x) :- 



1 



2m(2Z + 1) 2 2m(2Z + 1) + 



e(x) 



XB e (*) + 



dx - (2Z + l)biV, 



d R (x) 



where we have here also used a standard method to handle the non-square- 
integrability of the potential at infinity, by adding and subtracting a con- 
stant b > 0. Applying the Lieb-Thirring inequality (|59p with the resulting 
potential V% given above, we find 

(iP,H(N,R)Tj>) 

> -2mC' A v 2 (2Z + l) 2 J 2 (2m(2Z + 1) + a' 1 ) X B e + " 6 ) 

-(2Z + l)Mr4^1 
fc=i 



> -2mCA ^ 2 (2Z + I) 2 E ( / 4 (2m(2Z + 1) + 

fc=l \JB £k (R k ) 



£ k 1 f dx 



+ 







/ 


f, \,, »1 







da; 



Z 2 /(2Z + 1) : 
16mC' A v 2 D k 



{2Z + l)bN 
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K 



> -2mC' A v 2 {2Z + l) 2 ^ ( 4 {2m(2Z + 1) + e k v )\e 

1 



+ 2vr 



fe=l 

(r _1 — 6) 2 r dr 



UAmC' k v 2 D k 



(2Z + 1)6JV 



A 



> -4vrmCi i/ 2 (2Z + l) 2 £ ( 4 + l 6 ^|±^( £fc 6) 2 



+ 



fe=i 

1 3 



£fe6<1 ~ 288^^ (D * } 



(2Z + 1)WV. 



Finally, taking et '■= min{6 , D^} and in the latter case using that 

— Inx < [— lnc— 1]_|_ 

x 

for < x < 1 and c > 0, we have the further bound 



> -4irmC' A v l (2Z + l) 2 i^T 5 + 16 



,m 2 (2Z + If 



+ 



In- 



2887rmC" A 



□ 



- (2Z + l)bN 

> -Cv 2 m{2Z + l) 2 (K + JV) 

for some C > 0, where in the last step we chose b := v 2 m(2Z + 1). 

6.3 External potentials in the Heisenberg case 

Let us end with also considering an application of the above local exclusion 
and Lieb-Thirring inequalities for one dimensional statistics in the Heisen- 
berg case. We think of N identical such particles placed in an external con- 
fining potential V, and would like to obtain a lower bound for the ground 
state energy Eq. We have from Theorem 1181 that for a > 1 

>inf(/ {Cutu{a) 2 ~p{xf + V(x)p{x))dx : p : M -> M> , / p = ivl . 



Because of the local approximation p to the density, this minimization prob- 
lem is less tractable than e.g. the anyonic case in Section I6TT1 However, what 
we can do instead is to consider a partition V of the real line into a collection 
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of finite subintervals {Ij}jL 1 and an exterior domain I ext , and on each such 
interval apply Lemma [12] or the local form of Theorem [THl 



T H + / Vp > V f 2$ + f Vp]+ f Vp 




3=1 

where we define K- := inf r 1/ and 

j j j 

f(^) := max{0, pj - 1, X{ Pj >2} Pj/32}, Pj := / p. 

,//,■ 

Here we used (JMD and (|5SJ) with C H > 1/32, and note that p,- - 1 < p?/32 
and > 2 implies pj > p c , with the critical point p c 5.068 > 2. We have 
hence reduced to a finite-dimensional minimization problem 

E > E[V] := inf ^£ v [p] : p = p ext ) E R^ +1 , £ p, + p ext = AT 

:= X] ( fjp + ^ ) + ^xtPext, 

j= i V I il / 

for which a minimizer p exists due to the convexity of 8-p. Depending on N 
and the details of V, a hopefully good lower bound for Eq (and an estimate 
for the ground state density p) can then be found by maximizing over all 
suitable partitions V, 

E > sup{E[V] : V = ((IjO^iext) partition of R, M G N}. 

As an illustration of this procedure we can consider the following example, 
which in the harmonic case V(x) = ^|x| 2 can be compared to the exact 
ground state energy for the corresponding Calogero-Sutherland model [H 
[55] . Eq = |o;iV(l + a(N — 1)). We can also compare to the approximative 
Thomas- Fermi method discussed in |50[ 151]. 

Proposition 23. Given an external confining potential V(x) = c^\x\^ , with 
c, p > ; and assuming N S> (£ft(a)c) 2 /^ and ^u(a)cN 3> 1, we have for 
the ground state energy of H = Tr + Ylj=i V{xj) 

E > C(p,N)(£a(a)c)l** , (84) 



where C(p,N) > is a constant s.t. 

v r, m ( ^ + 2 nj) / r(|+j) 



H+2 
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In particular, liminfAr^oo Eq/N 2 > ^£h(oO w for p = 2 and c = u/y/2. 
Proof. Consider a simple partition of the form 

V = ([-{k + l)a, -fco]^ 1 , [ka, (k + l)a]%£, (-00, -Ma] U [Ma, +00)) , 

a > 0, for which we have (making use of the symmetry, and with a suitable 
relabeling of pj) 



M-l / p ( \2 \ 
£<p[p] = 2 J2[ ^T-C(Pk) + ^(kaypH + c^May Pext . 

1 — n V ^ / 



fc=0 

We extremize the functional 



F[p, A] := £ v [p] - A \J2 Pj + Pcxt + AiV 



= 2 E ( ^1r-£(pk) + (c"a"fc" - A) p k j + (cVM' 4 - A) Pcxt + AiV 
fc=o ^ a ' 

to find, if M is chosen large enough, p cxt = and pk = for c^a^k^ — A > 0, 
= 1 for < A - c^a^k^ < £,u(a) 2 /a 2 , and otherwise 



) aW32/3 r- — 

p k = max < p c , v VA - c^gW 

£hO) 



Hence, A r = 2 Pfc + Pcxt 



V 2 max I p c , Q y^ 2 ( 3 \/A - c^fc^ 

^ X 61(a) 



k>0 s.t. 
A- C M a MfeM>^ H ( a )2/ a 2 



E 



fc>0 s.t. 
0<A-c' i a' i fc^<$H(a) 2 /a 2 



Now, let a := 1/c and assume that N > (£c) 2 ^ and £ciV > 1. Then 
£cN ~ ^ 2 v / 32/3v / A- P 



fc>0 s.t. 
A-fc^V 



where 



/•A"" 

~ / VA - P dfc = 8^/273 /(p)A 1/2+1/M , 

jo 

/(/*) := /'vT 
J 



x' 1 (fx 



2 r(| + J) 
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( s/3/2 \ f+ 2 
Hence, A ~ I gjA ^cN 1 and 



E > £ v [p. 



fc>0 s.t. 



8x72/3 ^^i(A-F)f + F^j ~ 8^2/3 J(/.) A " 



3,1 



where 



2 2^ r(§ + i 



J( jU ) : =i /" (1 - a^)f cfe + / xVl - x^dx 
3 Jo Jo 

It follows under the above conditions that 

which is the asymptotics of the r.h.s. of (|84p . □ 
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